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Abstract 

Wc investigate strong coupling effects in a covariant massive gravity model, which is a candidate for a 
ghost-free non-linear completion of Fierz-Pauli. We analyse the conditions to recover general relativity 
via the Vainshtein mechanism in the weak held limit, and hnd three main cases depending on the 
choice of parameters. In the hrst case, the potential is such that all non-linearities disappear and the 
vDVZ discontinuity cannot be avoided. In the second case, the Vainshtein mechanism allows to recover 
General Relativity within a macroscopic radius from a source. In the last case, the strong coupling of 
the scalar graviton completely shields the massless graviton, and weakens gravity when approaching 
the source. In the second part of the paper, we explore new exact vacuum solutions, that asymptote to 
de Sitter or anti de Sitter space depending on the choice of parameters. The curvature of the space is 
proportional to the mass of the graviton, thus providing a cosmological background which may explain 
the present day acceleration in terms of the graviton mass. Moreover, by expressing the potential for 
non-linear massive gravity in a convenient form, we also suggest possible connections with a higher 
dimensional framework. 



1 Introduction 

Attempts to build a theory of massive gravity date back to the work by Fierz and Pauli (FP) in 1939 
[1]. They considered a mass term for linear gravitational perturbations, which is uniquely determined by 
requiring the absence of ghost degrees of freedom. The mass term breaks the gauge invariance of General 
Relativity (GR), leading to a graviton with five degrees of freedom instead of the two found in GR. There 
have been intensive studies in to what happens beyond the linearised theory of FP. In 1972, Boulware 
and Deser found a scalar ghost mode at the non-linear level, the so called sixth degree of freedom in the 
FP theory [2]. This issue has been re-examined using an effective field theory approach [3], where gauge 
invariance is restored by introducing Stiickelberg fields. In this language, the Stiickelberg fields acquire 
non-linear interactions containing more than two time derivatives, signalling the existence of a ghost. In 
order to construct a consistent theory, non-linear terms should be added to the FP model, which are 
tuned so that they remove the ghost order by order in perturbation theory. 

Interestingly, this approach sheds light on another famous problem with FP massive gravity; due to 
contributions of the scalar degree of freedom, solutions in the FP model do not continuously connect 



to solutions in GR, even in the limit of zero graviton mass. This is known as the van Dam, Veltman, 
and Zakharov (vDVZ) discontinuity [4, 5]. Observations such as light bending in the solar system would 
exclude the FP theory, no matter how small the graviton mass is. In 1972, Vainshtein [6] proposed a 
mechanism to avoid this conclusion; in the small mass limit, the scalar degree of freedom becomes strongly 
coupled and the linearised FP theory is no longer reliable. In this regime, higher order interactions, which 
are introduced to remove the ghost degree of freedom, should shield the scalar interaction and recover 
GR on sufficiently small scales. 

Until recently, it was thought to be impossible to construct a ghost-free theory for massive gravity that 
is compatible with current observations [7, 8] . A breakthrough came with a 5D braneworld model known 
as Dvali-Gabadadze-Porrati (DGP) model [9]. In this model there appears a continuous tower of massive 
gravitons from a four dimensional perspective, and GR can be recovered for a given range of scales, due 
to strong coupling interactions [10, 11]. In this paper, we explore the consequences of a promising new 
development along these lines that seem able to provide a consistent theory of massive gravity directly 
in four dimensions. In order to avoid the presence of a ghost, interactions have to be chosen in such 
a way that the equations of motion for the scalar degrees of freedom contain no more than two time 
derivatives. Recently, it was shown that there is a finite number of derivative interactions that give rise 
to second order differential equations. These are dubbed Galileon terms because of a symmetry under 
a constant shift of the scalar field derivative [12]. Therefore, one expects that any consistent non-linear 
completion of FP contains these Galileon terms in the limit in which the scalar mode decouples from the 
tensor modes, the so-called decoupling limit. This turns out to be a powerful criterion for building higher 
order interactions with the desired properties. Indeed, following this route, de Rham and Gabadadze 
constructed a family of ghost-free extensions to the FP theory, which reduce to the Galileon terms in the 
decoupling limit [13]. 

In this work, we investigate the consequences of strong coupling effects in this theory. We first re- 
express the potential for the most general version of non-linear massive gravity, as developed by de Rham, 
Gabadadze and collaborators, in a particularly compact, easy-to-handle form. Among other things, this 
way of writing the potential suggests intriguing relations with a higher dimensional set-up, which might 
offer new perspectives for analysing this theory. Moreover, we show that, for certain parameter choices, 
the potential for massive gravity coincides with the action describing non-perturbative brane objects, 
independently supporting connections with a higher dimensional framework. 

Armed with these tools, we then focus on the Vainshtein mechanism for this potential. We show 
that this theory is able to reproduce the behaviour of linearised solutions in General Relativity below 
the Vainshtein radius, but only in a specific region of parameter space. This result provides stringent 
constraints on non-linear massive gravity. Moreover, we are able to physically re-interpret these findings 
in the decoupling limit in terms of an effective theory with Galileon interactions. We show that the 
condition to successfully implement the Vainshtein mechanism is associated with the sign of a direct 
coupling between the massless graviton and the scalar degrees of freedom. 

We also present new exact solutions in the vacuum that asymptote de Sitter or anti de Sitter space 
depending on the choice of the parameters. Asymptotically de Sitter configurations can be expressed in 
an explicit time-dependent form. These solution may provide an interesting background for the observed 
Universe where the rate of the accelerated expansion of the Universe is set by the graviton mass. A 
small graviton mass, as required by the solar system constraints on deviations from standard General 
Relativity, is then in agreement with the observed value of the cosmological constant. On the other 
hand, asymptotically anti de Sitter configurations may have interesting applications to the AdS/CFT 
correspondence. 
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The paper is organized as follows: in Section 2, we discuss how to construct a non-linear potential for 
massive gravity and point out new connections with a higher dimensional set-up. In Section 3, we show 
how linearised Einstein's gravity is recovered within a certain macroscopic radius from a mass source, 
via the Vainshtein mechanism. Only a subset of parameter space presents a successful Vainshtein effect 
in the weak field limit. For a better understanding of the theory and in particular of the ghost mode, it 
is important to find analytic non-linear solutions; we present new exact solutions in Section 4. Finally, 
we conclude in Section 5, leaving technical details of the calculations developed in the main text to the 
Appendixes. 

2 Covariant non-linear massive gravity 

We start by introducing the covariant Fierz-Pauli mass term in four-dimensional spacetime 

C FP = m 2 ^gU {2 \ = (H^H^ - H 2 ) , (2.1) 

where the tensor H^ u is a covariantisation of the metric perturbations, namely 

9p,v = iJtMv + \v = H fw + d^du^rjap. (2.2) 

The Stiickelberg fields cj) a = {x a — tt 01 ) are introduced to restore reparametrisation invariance, hence 
transforming as scalars [3]. The internal metric r\ a p corresponds to a non-dynamical reference metric, 
usually assumed to be Minkowski space-time. Therefore, around flat space, we can rewrite H^ u as 

H^v = + rip^ir 13 + rj ail d u ir a - ^d^d^ ', 

= ( 2 - 3 ) 

From now on, indices are raised/lowered with the dynamical metric g^ v , unless otherwise stated. For 
example, Hi, = g^ p H pu . Moreover, the Lagrangian (2.1) is invariant under coordinate transformations 
x p — > x p + provided 7r^ transforms as 

7T^ -> 7T^ + (2.4) 

The scalar component tt of the Stiickelberg field can be extracted from the relation tt p = rj^ u d^ir/As, 
with A| = m 2 M p i (the meaning of the scale A3 will be explained in the following). The dynamics of 7r 
are the origin of the two problems discussed in the introduction: the DB ghost excitation and the vDVZ 
discontinuity. With respect to the first problem, as noticed by Fierz and Pauli, one can remove the ghost 
excitation, to linear order in perturbations, by choosing the quadratic structure h pu h pu — h 2 . When 
expressed in the Stiickelberg field language, by means of the scalar-graviton tt, higher derivative terms 
in the action are arranged in a such way to form a total derivative, leading to second order equations of 
motion. However, when going beyond linear order, the equation of n acquires higher time derivatives, 
signalling the presence of a ghost mode [3]. Remarkably, de Rham and Gabadadze were able to construct 
a potential, tuned at each order in perturbations, to give a total derivative for the dangerous terms, 
leading to equations of motion that are at most second order in time derivatives [13]. 

We now review their construction, introducing alternative ways to express the potential, which provide 
a new connection with a five dimensional point of view. In terms of the helicity zero mode, corresponding 
to the field ir, we can write the tensor H^ u of Eq. (2.3) as 

2 1 2 
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where = d^d^TT and n 2 ^ = 11^11". At a given order n in perturbations, the idea is to add terms of 
the form 



nr^-gU™ = m 2 v^ £ c?(^)(iT), (2-6) 



8=0 

to the FP action (2.1), and to choose the coefficients c n in order to get a total derivative for the leading 
contributions of the scalar mode, namely (ILj„) n . The key finding of [13] is that these total derivatives 
are unique at each order, and that the series stops at quintic order in perturbations. Let us describe in 
more detail the structure and origin of these terms. Following the notation of [13], the total derivatives 
are given by 

42 = - E(-ir^^(trn-)ct m \ (2.7) 

with C^} r = 1 and 

— t r TT 

42- = (trn MI/ ) 2 -trn^, 

4?r = (trn^) 3 -3(trn^)(trn^) + 2trn^, 

4tt = (trIV) 4 - 6(tr IV) 2 (tr IT 2 ,) + 8(tr n^)(tr II 3 v ) + 3(tr II 2 A 2 - 6trnJ„ . 

and 4er vanishes for n > 4. These expressions are related to a matrix determinant (see also Note Added 
at the end of this paper). To see this, consider a generic squared real matrix A, and a complex number 
z. Then, the following formula holds 



det(I + ,zA) = l + YVdet(A) (2. 



n=l 

where det n (A) can be written in terms of traces as 

det(A) = trA, 

det(A) = I ((trA) 2 -trA 2 

det(A) = -f(trA) 3 -3(trA)(trA 2 ) + 2trA 3 
3 6 \ 

det(A) = ^((trA) 4 -6(trA) 2 (trA 2 ) + 8(trA)(trA 3 ) + 3(trA 2 ) 2 - Glr.l 1 ). (2.! 



Moreover, all terms det n (A) with n > 4 vanish for a 4 x 4 matrix. Therefore, for the choice AJ^ = 11^ we 

get the simple relation 
determinants of the form 



get the simple relation = n\ det n (II), and the series indeed stops at n = 4. If one chooses a sum of 



^det(I + z 4 n) -4, (2.10) 



one can generate each det n (n) term with a separate coefficient (3 n , provided a solution to Yli=i z ? = Pn 
exists, which is guaranteed by the Newton identities. Then, the Lagrangian for the helicity zero mode 
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(that is, neglecting for the moment the contributions of tensor modes, and of the vector components of 
the Stuckelberg field) is 

4 

Ar = y>„det(II). (2.11) 

n=l 

We now briefly turn away from the present discussion, and show an interesting way to construct 
the Lagrangian (2.11) from a higher dimensional point of view. Consider a five dimensional Minkowski 
spacetime, and embed on it a test 3-brane (i.e. we do not include back-reaction from brane dynamics). 
Under this assumption, the five dimensional Riemann tensor vanishes; using the Gauss equation, the 
intrinsic curvature on the brane is related to the extrinsic curvature as (R a p^$ is constructed in terms of 
brane induced metric) 

R a MS = K«Kp S -KfK M , R ilv = KK ilv -KlK vay R = K 2 - K^K^ . (2.12) 
We then consider a four dimensional Lagrangian given by 

£brane = V~9 a l K + + a^GB + OL 4 R GB , (2.13) 

where R is the Ricci scalar, Rgb is the Gauss-Bonnet term, Rqb = R 2 — ^R 2 ^ + R^vapi K is trace of 
the extrinsic curvature, and the /Cqb is the boundary term associated with the five dimensional Gauss- 
Bonnet term: fC GB = K 3 - 3KK 2 U + 2K 3 U . Using the expression for the intrinsic curvature in terms of 
the extrinsic curvature, Eq. (2.12), the Lagrangian (2.13) can be written as 

4 

Arane = " V/3ndet(JC), (2.14) 
* — • n 

n=\ 

where j3 n = —n\a n , thus it has exactly the same structure as Eq. (2.11). We then denote the position 
modulus of the probe 3-brane as tt. The induced metric on the brane is determined by n as 

9ij,u = riiLv + d^Trd u TT, (2.15) 

and the extrinsic curvature is given by 

^ = 7^vr, 1= —=L==. (2.16) 
ya + (dvr) 2 

If we take a limit Ott <C 1, i.e. 7 — > 1, we find that the extrinsic curvature is simply K^ u = Then the 
Lagrangian (2.14) reduces to (2.11). This suggests that there may be a higher-dimensional interpretation 
behind the Lagrangian (2.11). Although we find these arguments very compelling, so far we have not 
been able to pursue these connections further, and for this reason we will not develop them in this work. 
But we should note that what we discussed follows the same construction as in the so called DBI Galileon 
[27]. If we do not take the limit 7 — > 1, the Lagrangian (2.14) becomes non-trivial, and it reproduces 
the Galileon terms. Since the four dimensional Gauss-Bonnet piece is a total derivative, there is no 
contribution from this term even away from the 7—7-1 limit. 

After this digression to five dimensions, we would like to return to Lagrangian (2.11), and discuss 
how to render it fully covariant. To do so, we need to understand how to go back from the field tt to the 
original Stuckelberg fields 4>^, in order to restore the dependence of the vector mode and the full metric. 
Here we follow the approach discussed in Ref [14] . If only the scalar mode ir is considered, then we can 
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solve for Ti^ u in terms of (f)^, using (2.2) and (2.5). The result is a second order algebraic equation for 



IV, with solution iy = A 3 5^ - (y {d<j>)r]{d<l>) T j , where [{d^dcpf}^ = d^ a d v ^r] a p, and just 
to remind the reader A| = m 2 M p i. The previous expression is written in terms of the square root of a 



matrix, formally understood as \M U \M a = M". 



We now go beyond the pure scalar sector case, and define 

V = s; - (y/WMdt)^ " = 5; - (Vi-s- 1 ^ = - (Vg-'lv + Q])^ ■ (2.17) 

Notice that the previous quantity contains also contributions from vector and tensor degrees of freedom. 
On the other hand, by construction fCu becomes liu /A3 when only the scalar mode is considered. The full 
non- linear Lagrangian for massive gravity is then constructed by substituting II by K, in (2.11). Namely 



a x det(/C) + 2a 2 det(/C) + 6a 3 det(/C) + 24a 4 det(/C) (2.18) 



where a n = —n\j3 n , and the determinants det n (/C) are defined in Eqs. (2.9) with A = K.. The second 
term, for positive 02, reduces to the Fierz-Pauli term (2.1) when expanding the Lagrangian in terms of 
H^ u around the Minkowski metric. Therefore, since we would like to have the Fierz-Pauli as the first 
correction to Einstein's gravity at leading order, we are not be interested on the contributions from the 
first term, deti tC. Then, from now on we set ct% = 0. As a result, a family of non-linear massive gravity 
Lagrangians can be written as 

M 2 

£ = — SL (R - 2A - m 2 U) , (2.19) 



2 

where IA = Cfc with ct\ = and a% = 1. These Lagrangians are parametrised by m, 03 and 04; moreover 
we added a bare cosmological constant A. 

In this family of Lagrangians, there is a special choice of parameters: 03 = —1/3 and 04 = 1/12. It 
corresponds to the choice z = — 1 in the expansion (2.8) with A = K,. The Lagrangian is 



C NG = 2m 2 7^(det(I- K) - tr/C) = 2m 2 (yj- detid^d^r]^) - tr K\ , (2.20) 

where the first term is the Nambu-Goto type of action for a bosonic 3-brane [15]. However, we will discuss 
in what follows that this particular choice, even though has a striking physical interpretation, does not 
allow to recover the GR solutions via the Vainshtein mechanism. 

The «2 term in the the potential (2.18) was first suggested in [14] as non-linear completion of FP 
theory. It was shown that this term reduces to a particular choice of Galileon terms in the decoupling 
limit, without any coupling between the scalar mode and the massless graviton. Formally, the decoupling 
limit corresponds to a limit in which the scale A3 = m 2 M p i is kept fixed, while sending M p i — > 00 and 
m — > 0. Once 03 and 04 are included, the picture changes in an interesting way, and couplings between 
7r and the massless graviton appear, even in the decoupling limit. These couplings have important 
theoretical and observational consequences; as we will discuss in detail in the next Section. Interestingly, 
these mixing terms are finite in number, and do not spoil the fundamental property that the equations of 
motion for ir are second order [13]. It is still an open question whether this remains true away from the 
decoupling limit, ensuring the absence of ghost degrees of freedom. A full Hamiltonian analysis should 
be carried out to set a final word on the subject. However, there are hints that the theory is ghost-free 
perturbatively, and for the particular choice of 03 = = this has been shown up to and including 
quartic order in perturbations [14]. 
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We have now the necessary ingredients to discuss the second problem addressed in the Introduction, 
namely the vDVZ discontinuity. 



3 Vainshtein mechanism at work 

In [16], we showed that Vainshtein mechanism applies for Lagrangians as (2.19), setting to zero the coef- 
ficients 0:3 and o;4, and the bare cosmological constant. (See Ref. [18] for spherical symmetric solutions in 
the FP theory). Here we extend the analysis to arbitrary coefficients. We determine stringent constraints 
on the parameter space of these theories, in order to ensure that the Vainshtein mechanism works. As 
we are going to discuss, our results find a natural interpretation in terms of the dynamics of helicity-zero 
mode, in the decoupling limit. 

In order to discuss solutions associated with Lagrangian (2.19), it is convenient to express /C, given 
in Eq. (2.17), in terms of matrices as 

K = I-VM, (3.1) 

where I denotes the identity matrix and M. = g^ 1 [ry + Q]. The task is to calculate the trace of A4 n . 
Given that Ai is a square matrix, it is diagonalizable, and can be expressed as Ai = UDIA^ 1 , for some 
invertible matrix U, where D is a diagonal matrix containing the eigenvalues of A4. We shall call the 
eigenvalues Ai, . . . , A4. Then, since Ai n = U D n U , the traces in the formulae above can be expressed 
in terms of eigenvalues 

tvM n = Y,K, (3-2) 

i 

and the traces of K, n result 

tr/C = 4-trVM 

tr/C 2 = 4 - 2tr V / A4 + trA4, 

tr/C 3 = 4-3tr\/A? + 3trAf -trX 3/2 , 

tr/C 4 = 4-4trv / A4 + 6trA4 - 4tr A4 3/2 + tr A* 2 . (3.3) 

Using the formulae for expressing the determinants in terms of traces (2.9), we can easily construct the 
potential. 

We now discuss the conditions to recover GR results in the small graviton mass limit, within a certain 
radius from a mass source. In particular, we are interested to determine under which circumstances the 
Vainshtein mechanism applies. In order to do so, we study spherically symmetric perturbations around 
flat space, expressed in spherical coordinates as ds 2 = — dt 2 + dr 2 + r 2 d£l 2 , with dQ 2 = d6 2 + sin 2 Odcf 2 . 
We start our discussion using the unitary gauge, 7r M = — ^ = 0. Consider the following Ansatz for 
the metric 

ds 2 = -N(r) 2 dt 2 + F{ry x dr 2 + r 2 H ( r y 2 dfl 2 , (3.4) 
that reduces the potential in (2.19) to 

r 2 sin((9) 



-gU 



FH 2 



1 2 VF((2H -3)N + 1) + H 2 N + H(2-6N) + 6N -3 

-6a 3 (H - 1) \y/F((H - 3)N + 2) - 2HN + H + AN — 3 
-24q 4 (1- Vf)(1- H) 2 (l- AT) J. (3.5) 
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Notice that in GR one can set H(r) = 1 by a coordinate transformation, but this is not possible here, 
since we have already chosen a gauge. The field equations are obtained by varying the action (2.19) 
with respect to N, H and F. The TV-equation is the Hamiltonian constraint, so it only depends on 
F and H. Since the equation for H is quite complicated, we instead consider a combination of the 
three equations which gives V^G^ = 0, where G^ u is the Einstein tensor; it corresponds to the Bianchi 
identity. Therefore, we work with the Hamiltonian constraint, the Bianchi identity and the equation for 
F. The corresponding expressions are lengthy, so we relegate them to Appendix A. 
First, let us study solutions in the weak field limit, by expanding N, F and H as 

N = l + n, F = l + f, H = l + h, (3.6) 

and truncating the field equations to first order in these perturbations. As we will see in what follows, 
this linearisation procedure is not completely consistent for all values of the radial coordinate r, and 
we will need to improve it. In order to analyse the system, it is convenient to introduce a new radial 
coordinate 

H (r) 

so that the linearised metric is expressed as 

ds 2 = -(1 + 2n)dt 2 + (1 - f)dp 2 + p 2 dn 2 , (3.8) 

where / = / — 2h — 2pb! and a prime denotes a derivative with respect to p. As discussed above, one 
should be careful with this change of coordinates, since, after fixing a gauge, a change of frame in the 
metric modifies the Stiickelberg field ir^ as well. However, for the moment, let us focus on the change of 
the metric part; later we will discuss what happens to tt^. At linear order, the equations for the functions 
n(p), f(p) and h(p) in the new variable p are 

= (m V + 2) / + 2p ff' + m 2 p 2 ti + 3m 2 phj , (3.9) 

= m 2 p 2 (n - 2h) - 2pn' - f, (3.10) 
= f + pn. (3.11) 

In this linear expansion, the solutions for n and / are 

2n = -^e"-, 
3p 

4GM 

f = -^—{l + mp)e- m ^ (3.12) 
3p 

where we fix the integration constant so that M is the mass of a point particle at the origin, and 
8vrG = M~ 2 . These solut ions exhibit the vDVZ discontinuity, since the post-Newtonian parameter 
7 = f /2n is 7 = 1/2(1 + mp), which in the massless limit reduces to 7 = 1/2, in disagreement with GR, 
and with Solar system observations (7 = 1 in GR, while observations provide 1 — 7 ~ 10~ 5 [17]). 

However, in order to understand what really happens in this limit, we must also analyse the behaviour 
of h as m — > 0. For doing this, we consider scales below the Compton wavelength mp <C 1, and at the 
same time ignore higher order terms in GM. Under these approximations, the equations of motion can 
still be truncated to linear order in / and n, but since h is not necessarily small, we have to keep all 
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non-linear terms in h. The resulting equations are then (see Appendix A for their derivation) 
= 2/ + 2pf + m 2 /? 2 ( [1 - 2(3q 3 + l)h + 3(a 3 + 4a 4 )/i 2 ] [(2 + f)pti + (1 + h)f 

+6h[l - (3a 3 + l)h + (a 3 + 4a 4 )/i 2 ]) , (3.13) 

= -/ - 2pra' + m 2 p 2 (n-2[l + n+ (3a 3 + l)n]h + [(3a 3 + l)(n + 1) + 3(a 3 + 4a 4 )n]/i 2 ) , (3.14) 

= pn' [-1 + 2(3a 3 + l)h - 3(a 3 + 4a 4 )/i 2 ] - f[l - (3a 3 + l)h] . (3.15) 

We start with the iV-equation (3.13). Since it only depends on / and h, one can solve for / in terms 
of h, including all non-linear terms in h. The solution, dropping higher order terms proportional to m 4 , 
{GMf and m 2 GM, is 

/ = -2^- - (mp) 2 \h - (1 + 3a 3 )h 2 + (a 3 + 4a 4 )/i 3 l . (3.16) 
P 

Then we take the second equation (3.14), obtained by varying the action with respect to F, and use the 
solution (3.16) for /. We find an expression for n as a non-linear function of h. It turns out to be simpler 
to work with n', given by 

2C M 

2pn' = (mp) 2 \h - (a 3 + 4a 4 )/i 3 l , (3.17) 

P 

where again we have dropped terms with m 4 , (GM) and m 2 GM. Finally, the constraint equation (3.15) 
gives an equation for h, after substituting the solutions for n and / given respectively by (3.17) and 
(3.16). We should stress that this equation for h is exact, so there are no higher order corrections. It is 
given by 

(rlf [1 - 3(a 3 + 4a 4 )/i 2 ] = -(mpf^h - 3(1 + 3a 3 )h 2 + [(1 + 3a 3 ) 2 + 2(a 3 + 4a 4 )] h 3 

--(ag + 4a 4 ) 2 /i 5 }- (3.18) 



P 



If we linearise the equations (3.16), (3.17) and (3.18) with respect to h, we recover the solution (3.12) 
on the other hand, below the so-called Vainshtein radius 



PV - (^Y 3 , (3.19) 



\ m z J 

h becomes larger than one. Therefore, in this regime, we have to include higher order contributions due 
to h to equations (3.16) and (3.17). There are three qualitatively different cases, depending on the values 
of the parameters a 3 and a 4 : 

• Case with a 3 = — 1/3, a 4 = 1/12. This is a special situation, since all higher order contributions 
in h vanish from equations (3.16)-(3.18). Therefore, there is no Vainshtein effect and the solutions 
to the equations are those given in (3.12) for p < 1/m. This model, corresponding to the bosonic 
3-brane Lagrangian (2.20), is ruled out by Solar system observations. 

• Case with a 3 = — 4a 4 ^ —1/3. For p <C pv we can solve for h from the last equation (3.18) keeping 
only the highest order terms in h. Then, the solution for h is 

h = - ( i-L t ^ T- (3 - 20) 
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which implies \h\ 3> 1 for p <C py, as expected. We can then use this solution and equations (3.16) 
and (3.17) to give the expressions for n and / within the Vainshtein radius, namely 



2n 



f 



2GM 
P 

2GM 
P 



1 



1 



1 



2(1 - 12a 4 ) 2 /3 \pv J 



2(1 - 12a 4 ) 1 /3 \ pv J 



(3.21) 



Therefore, the corrections to GR solutions are indeed small for p smaller than the Vainshtein radius 
PV- 

Case with ^ — 4a 4 . This case can be divided in two: 03 + 4a 4 > and Q3 + 4a 4 < 0. The latter 
is the most intriguing case, and we focus on it. In the limit p <C py, the solution for h is given by 



«3 + 4a 4 



1/3 



PV 
P 



2(l + 3a 3 ) +3(a 3 + 4a 4 ) 



2 (a 3 + 4a 4 ) E 



_P_ 

PV 



(3.22) 



so that \h\ 3> 1. Notice that, in solving Eq. (3.18), we include also the next-to-leading order, that 
results to be linear in p. It turns out that this expression for h provides a correction of the same 
order as GM/p in the n and / equations. Indeed, plugging this expression in the equations for n 
and /, one gets 



2n 
f 



O 



p \2\ GM 
Pv' 



P 



p\GM 

W' p 



(3.23) 



Surprisingly, the contribution from the scalar mode h exactly cancels the usual 1/p potential at 
leading order. So, gravity becomes weaker approaching the source, for distances smaller than 
Vainshtein radius, and larger than Schwarzschild radius. This implies that the strong coupling of 
the scalar graviton not only shields interactions of the scalar mode h, but also those of the massless 
graviton. As we will see later while analysing the decoupling limit, this is related to a coupling 
between the scalar mode and graviton that cannot be removed by a local field transformation if 
03 7^ — 4a 4 [13]. In this case, local tests of gravity would completely rule out the predictions of the 
theory at leading order in h^ v . However, when 03 + 4a 4 > 0, apart from the previous solution, a 
second solution exists in the limit p <C pv [25], where h = 1/ \J 3(«3 + 4a 4 ) + 0(p 3 ). In this case, 
the metric fields n and / present small corrections to the GR solution within the Vainshtein radius, 
and the solution cannot be ruled out by Solar system tests. 

In order to fully understand the system with 03 7^ — 4a 4 , it is imperative to study how the two 
possible solutions extend beyond the Vainstein radius, both in vacuum and in the presence of 
matter. This is beyond the scope of the present analysis and will be studied in a future publication 
[26]. 

To summarize, only the choices ce3 + 4a 4 > and = — 4a 4 7^ —1/3 allow to recover standard GR in the 
weak field limit, below the Vainshtein radius py . This fact imposes a stringent constraint in parameter 
space for the theory described by Lagrangian (2.19). 
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Now, we would like to go back to the issue of the coordinate transformation introduced in Eq. (3.7), and 
discuss another way to interpret the previous results. As mentioned earlier, a coordinate transformation 
introduces a change in tt of the form (2.4). When changing r/H = p, we excite the radial component 
of the Stiickelberg field, as tt p = —ph. Thus the strong coupling nature of h is encoded in tt p = 
r\ pp dpTr/A^, when working with the coordinate p. As a result, the non- linear analysis previously done 
is more transparent in the decoupling limit, as developed in [13]. As previously mentioned, this limit is 
achieved by taking m — > and M p i — > oo, while keeping A3 = (m 2 M p ;) 1//3 fixed. This implies that, when 
substituting back into the full Lagrangian (2.19), we do not consider the vector mode, and expand 
the potential at linear order in haw The resulting expression describes the theory in the decoupling 
limit, and contains the quadratic Hilbert-Einstein piece, total derivatives in tt given by (2.11), and finally 
mixing terms between hav and tt. The mixing terms are 



with 



Yx« 



) 

fJ,U 



dh pv 



(3.24) 



and each J" is of order 0(IT). There are only three mixing contributions, which can be absorbed 
(except for X^) in the remaining terms by the non-linear field redefinition 



hnv = + ■zrz-f]av ; 5 ; 7 dpTrd u TT. 

1V1, - 



pi 



AlM pl l " 



After this field redefinition, the resulting Lagrangian is ([IT™] = trlP) [13] 



(3.25) 



C 



3(l + 3a 3 ) 



2A| 



{dix) 2 nir + 



2A| 



:i + 3a 3 ) 2 + 2(a 3 + 4a 4 ) ([n 2 ] - [n] 2 ) 



~(1 + 3a 3 )(a 3 + 4a 4 )(avr) 2 ([n] 3 - 3[H][n 2 ] + 2[H 3 ]) + ¥*hTX® 



PuV 1 



(3.26) 



where jCgr is the quadratic Einstein-Hilbert action for hu U , and 



uu 



-i(a 3 + 4a 4 ){6n 3 , - 6(tr n)n 2 , + 3 [(trn) 2 - (tr n 2 )] 

- [(trn) 3 - 3(trn 2 )(trn) + 2(trn 3 )] 7^} . 



(3.27) 



Interestingly, the Lagrangian obtained so far in the decoupling limit allows us to re-interpret the con- 
straints on a 3 and q 4 obtained earlier in this Section. One can observe that the remaining direct coupling 
between the scalar tt and graviton hav either vanishes or is negative for the two cases which successfully 
implement Vainshtein mechanism. 

Observe also that the kinetic terms for tt in (3.26) are precisely the Galileon terms, which give rise 
to second order differential equations for tt [12]. The non-linear structure of these terms involving tt 
is essential to recover GR in some range of scales distances; for the choice a 3 = — 4a 4 = —1/3, we 
obtain an Einstein frame Lagrangian for Brans-Dicke gravity with a vanishing Brans-Dicke parameter 
and no potential. However, this particular choice is not compatible with observations, as we have already 
mentioned during our general analysis, in the first part of this Section. 

We would like now to further analyse the equations of motion for this system in the decoupling limit. 
The Lagrangian (3.26) is exact in the decoupling limit, as there are no higher order terms besides those 
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shown, so the linearised Einstein equation for h^v is 



6G^(h)-j± XE X$=0. (3.28) 



(3) 

If we assume spherical symmetry, then it can shown that X^J is simply 

X $=-<*±^±^y, xg)=0. (3.29) 

Therefore, the solutions for the linearised metric, hu = —2h and h pp = —f, reduce to 

2GM p 2 /7r'\3 
p A^M p i V p) 

2pn = -/. (3.30) 
On the other hand, the equation of motion for tt derived from the action (3.26), is given by (see [12]) 

6 , , fn'V 2 r. v , ,1 /vr /x 3 



3(-J +^3(l + 3a 3 ) J^-J + ^| [(1 + 3q 3 )' + 2(q 3 + 4q 4 ) 
6Mpi(a3_+4a4) ^/V\ 2 _ M 



p 2 A| vr ; V P J ^P 3 M pl ' 

(3.31) 

where the integration constant is again chosen so that M is a mass of a particle at the origin. Using 
the relation between tt and h, h = —it'/ (m 2 M p ip), it is simple to check that the solutions for /, n and 
h given by Eqs. (3.16)-(3.18) agree with the expressions in Eqs. (3.25), (3.30) and (3.31). This confirms 
that the results obtained earlier in this Section are in perfect agreement with what is found in terms of 
the dynamics of the scalar field tt, in the decoupling limit. 

In summary, the Vainshtein mechanism applies for a 3 + 4a 4 > and a 3 = — 4a 4 ^ —1/3. Only for 
these choices the weak field GR results are fully recovered at distances smaller than the Vainshtein radius. 
If these solutions match the asymptotic configurations (3.12) beyond the Vainshtein radius, then there 
would be three phases: on the largest scales beyond Compton wavelength, m^ 1 -C p, the gravitational 
interactions decay exponentially due to the graviton mass, see Eq. (3.12). In the intermediate region 
pv < p < m" 1 , we obtain the 1/r gravitational potential but Newton's constant is rescaled, G — > AG/ 3. 
Moreover, the post-Newtonian parameter reduces to 7 = 1/2 in the m — > limit, instead of 7 = 1 of 
GR, showing the vDVZ discontinuity. Finally, below the Vainshtein radius p < py, the GR solution 
is recovered due to the strong coupling of the n mode (see, for example, Eq. (3.21) for the case a 3 = 
— 4a 4 7^ —1/3); in this regime the theory can then be rendered compatible with observations. 

The solution discussed here provides a testing arena for studying the Boulware-Deser ghost. Instead 
of expanding the action in H^ u around Minkowski spacetime up to higher orders in perturbations, we have 
the possibility to study linear perturbations around this non-perturbative solution, using the complete 
potential in Eq. (2.18). In order to obtain the full non-linear solution, matching the three phases we have 
described, a numerical approach is necessary. In the next section, we consider a different family of vacuum 
solutions for this theory, which can be obtained analytically, and can lead to interesting candidates for 
realistic backgrounds. 
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4 Exact solutions 



As we learned in the previous section, an essential property of this theory of massive gravity is the 
strong coupling phenomenon occurring in the proximity of a source. This allows, for certain regions of 
parameter space, to recover linearised General Relativity at sufficiently small distances by means of the 
Vainshtein mechanism. This behaviour, accompanied by the fact that Birkhoff theorem does not apply 
in this context, suggests that exact solutions for this theory, even imposing spherical symmetry, might 
be very different from the GR ones. 

In this section, we will exhibit new spherically symmetric exact solutions in the vacuum for massive 
gravity, that generalize the ones of [19] and [16]. In an appropriate gauge, the solutions are asymptotically 
de Sitter or Anti-de Sitter, depending on the choice of parameters. 

While in [16] we focused on the case a 3 = ce 4 = 0, we now generalize the analysis to arbitrary couplings 
in the Lagrangian (2.18). We adopt the unitary gauge and allow for arbitrary couplings 04, i = 2,.., 4 
(as explained earlier, we set to zero the coefficient a\). We start with the following form for the metric 
(for convenience, we implement slightly different conventions with respect to the previous section) 

ds 2 = -C(r) dt 2 + 2D(r) dtdr + A(r) dr 2 + B{r) dti 2 (4.1) 

so that, even though the spacetime is spherically symmetric, the metric contains a cross term dtdr. We 
choose the following Ansatz for the metric functions [19, 16], 



B(r) = b r 2 , 

C{r) = c + ^ + c 2 r 2 , 

A{r) + C{r) = Q , 
D 2 {r) + A{r)C{r) = A , (4.2) 

and use the equations of motion to fix the constant parameters &0j Q)> ci, C2, Qo, Ao- Einstein equations 
read 

G„ v = SttGT^ (4.3) 

with energy momentum tensor T^ u = — -^==|^r, and Ck given in Eq. (2.18). 

In General Relativity, diffeomorphism invariance allows one to choose the function B(r) to be B(r) = 
r 2 , so that bo = 1. In this theory of massive gravity, after having fixed the gauge, this choice is no longer 
possible and the equations of motion determine 60 • in order to do this, one observes that the metric Ansatz 
(4.1) leads to the following identity between components of the Einstein tensor: C{r)G rr + A{r)Gtt = 0. 
This combination on the energy momentum tensor provides the following value for 60, 



1 + 6a 3 + 12a 4 ± ^/l + 3a 3 + 9c^ - 12a 4 
3(1 + 3a 3 +4a 4 ) 



(4.4) 



The upper branch generalizes the result of [16], while the lower branch is specifically associated with 
theories in which a 3 and/or ct\ are non- vanishing. After Plugging the metric components (4.2) in the 
remaining Einstein equations, one can find the values for the other parameters. The corresponding general 
expressions are quite lengthy, and for this reason we relegate them to Appendix B. As a concrete, simple 
example, in the main text we work out the special case we focussed on the previous section, a 3 = — 4a 4 . 
In this solution is given by the following values for the parameters 
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bo 

co 

C-2 



4/1- 12a 4 \ 2 



9 V 1 - 8a 4 
Ao 

bo' 

m 2 A 
4(12a 4 - 1)' 



16(1 -12a 4 ) 4 + 81(1 -8a 4 ) 4 A 

CJo = 9 • (4.5) 

36[l + 4a 4 (-5 + 24a 4 )] 2 

The previous solution is valid for a 4 in the ranges a 4 < 1/12 and a 4 > 1/8. Notice that the case 

a 4 = 1/12 corresponds exactly to the Lagrangian (2.20), discussed in Section 2. We find that c\ and 
Ao are arbitrary; this vacuum solution is then characterized by two integration constants. The resulting 
metric coefficients can be rewritten in the following, easier-to-handle form: 

9 A /l-8a4\ 2 IM _ / ■. 4 /l-12a 4 \ 2 2 

A(r) = -A — — [p(r)+ 7 + l] , B(r) - 1 ] " 







- 8a 4 




12a 4 






- 8a 4 




12a 4 


9A / 


' 1- 


8a 4 
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8a 4 / 



C(r) = ~ A ( ; " n 7 ) [l-p(r)] , (4.6) 
with (// = — ci/co) 

p(r) = g + ^ " ^ 7= Jg- p-^ V-i (4 7) 

P ^- r + 9(1 -8a 4 ) 2 ' 7 ~ 81A V l"8a 4 7 1 j 

In order to have a consistent solution, we must demand that the argument of the square root appearing 
in the expression for D(r), Eq. (4.6), is positive. A sufficient condition to ensure this is that \x > 0, and 

r— 4 /l - 12a 4 \ 2 /A . 

< < - I — - 4 ) . (4.8) 



The metric then could be rewritten in a more transparent diagonal form, by means of a coordinate 
transformation. However, a coordinate transformation of the time coordinate is not permitted, since 
until this point we have adopted the unitary gauge. Therefore, we now renounce to this gauge choice, 
and allow for a non-zero vector tt^ of the form tt^ = (7To(r), 0, 0, 0). One finds that then the metric can 
be rewritten in a diagonal form, as 

ds 2 = -C(r)dt 2 + A(r)dr 2 + B{r) dtt 2 , (4.9) 

while the equations of motion for the fields involved are solved by 

KJ 9\l-8a 4 J l-p(r) ' 0W l-p(r) ' V ' 
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with C(r) and B(r) being the same as in Eq. (4.6). If one makes a time rescaling 



t 



4(1 - 12a 4 ) 2 
1 - 8o 4 ) 2 



9A 



1/2 



■f 



the resulting metric has then a manifestly de Sitter-Schwarzschild, or Anti-de Sitter-Schwarzschild form. 
This depends on whether 04 is smaller or larger than 1/12, as can be seen inspecting the function p(r) in 
Eq. (4.7). On the other hand, we should point out that this time-rescaling cannot be performed, without 
further introducing a time dependent contribution to ttq. As expected, the metric in Eq. (4.9) can also be 
obtained by making the following transformation of the time coordinate dt = dt + ir' dr to the original 
metric (4.1). This produces a non-zero time component for 7r^, that does not vanish even in the limit 
m — > 0. 

To summarize so far, we found vacuum solutions in this theory that are asymptotically de Sitter 
or Anti-de Sitter, depending on the choice of the parameters (Another family of solutions with similar 
behaviour, but obtained for different choices of parameters, have been recently discussed in [24]). Let us 
point out that it is also possible to include a bare cosmological constant term \J—g A to the Lagrangian 
(2.19). Our solutions to the Einstein equations, with our metric Ansatz, remain formally identical. The 
only difference is that the function p(r) in Eq. (4.7) becomes 



p (1 ~ 1204) 
PK ' r 9(l-8o 4 ) 2 



m 1 + - (1 - 12a<4) A 
3 



(4.11) 



Notice that the additional integration constant Ao can not be used to 'compensate' the contribution of 
the bare cosmological constant A via a self-tuning mechanism, since Ao does not explicitly appear in 
the previous formula. For asymptotically de Sitter solutions, 04 < 1/12, choosing \i = 0, the metric can 
also be written in a time dependent form, at the price of switching on additional components of 7r M , as 
mentioned earlier. After dubbing 



1 



m 



(1 - 12o 4 ) 



m 2 + - (1 - 120:4) A 



we can make the following coordinate transformation t = Ft (t, p) and r = Ft (r, p) with 



Ft(r,p) 
F r (r,p) 



3 Aq /2 m V 1 
3/ 1-80:4 \^ r / 2 



l-12o 4 \ , /sinh(^) + ^e"/ 2 
J arctanh 1 



\cosh(^) 



2 VI - 12o4< 

The metric becomes that of flat slicing of de Sitter 



ds 2 



where the Hubble parameter is given by 



2 



-dr 2 + e™ r (dp 2 + p 2 dn 2 ) , 



m 2 + - (1 - 120:4) A 



2(1 - 12o 4 ) 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



The Stiickelberg fields tx^ are now given by tt^ = (tt t (t,p), tt p (t,p), 0, 0), with tt t = ttq + r — F t , 
tt p = p — F r . Interestingly, the value of the Hubble parameter is ruled by the mass of the graviton; in the 
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case of vanishing bare cosmological constant, we have a self-accelerating solution, in which the smallness 
of the observed cosmological constant can be associated with the smallness of the graviton mass. 

This self-accelerating solution, which reduces to that found perturbatively in Ref. [20], appears as 
an ideal background to explain present-day acceleration. Notice that this configuration is remarkably 
similar to that in the DGP braneworld model [21], though there are important differences. In order to 
study the viability of our non-perturbative solution, it is necessary to study the behaviour of fluctuations, 
to confirm that there is no ghost. On the other hand, in the DGP model, the self-accelerating solution 
suffers from a ghost instability [10, 11, 22], which is related to the ghost in the FP theory on a de Sitter 
background. 

5 Discussion 

In this work, we investigated the consequences of strong coupling effects for a theory of non-linear 
massive gravity, developed by de Rham, Gabadadze and collaborators. We first re-expressed the complete 
potential for this theory in a particularly compact and easy-to-handle form. Among other things, this 
way of writing the potential suggested intriguing relations with a higher dimensional set-up, that might 
offer new perspectives. Moreover, we showed that, for certain parameter choices, the potential for massive 
gravity coincides with the action describing non-perturbative brane objects, independently supporting 
connections with a higher dimensional framework. 

We then studied the conditions to implement the Vainshtein mechanism in this context. The theory 
is able to reproduce the behaviour of linearised General Relativity below a certain scale, but only in a 
specific region of parameter space. This result provides stringent constraints on this non-linear massive 
gravity models. Moreover, we were able to physically re-interpret our findings in the decoupling limit, 
in terms of an effective theory with Galileon interactions for the scalar-graviton. We showed that the 
condition to successfully implement Vainshtein mechanism is associated with the sign of a direct coupling 
between the massless graviton and scalar degree of freedom which cannot be removed by a local field 
transformation. 

We also presented new exact solutions in the vacuum for this theory that asymptote either de Sitter or 
anti de Sitter space, depending on the choice of the parameters. Asymptotically de Sitter configurations 
can be expressed in an explicit time-dependent form, providing an interesting background for the observed 
Universe, where the rate of the accelerated expansion of the Universe is set by the graviton mass. A 
small graviton mass, as required by the solar system constraints on deviations from standard General 
Relativity, is then in agreement with the observed value of the cosmological constant. On the other hand, 
asymptotically anti de Sitter configurations may have interesting AdS/CFT applications. 

Our results naturally lead to various important questions to be future examined. It would be in- 
teresting to put in a firmer basis the connection between this theory of non-linear massive gravity and 
the higher dimensional set-up described in Section 2. Exploiting this relation might also shed light on 
the absence of ghost excitations at all orders in perturbations about Minkowski. In Section 3, we found 
the intriguing result that, in some region of parameter space, the coupling between the scalar degrees of 
freedom and graviton shields the interactions of linearised graviton so that gravity appears to become 
weaker as the source is approached. To understand in full detail the system in this case, it is however 
necessary to analyse the dynamics of higher order metric fluctuations and their couplings with the scalar 
sector. Although it seems unlikely that strong coupling dynamics of gravitational modes behave such to 
mimic General Relativity in relevant regimes (by producing a sort of Vainshtein effect at higher order 
in perturbations), we cannot exclude this possibility for all ranges of parameters. This paper provides 
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all the necessary equations to study this deep issue, away from the decoupling limit. It would also be 
interesting to study in more detail the exact vacuum configurations discussed in Section 4, in particular 
to understand whether de Sitter solution controlled by the graviton mass can be rendered stable under 
perturbations. If so, it can be considered as a serious candidate for the observed universe where the 
present day acceleration of the Universe is only due to gravitational degrees of freedom. 

Note Added 

While this work was in the last stages of preparation, two papers [23, 24] appeared, with some overlap 
with our Section 2 with respect to the formulation of the potential for non-linear massive gravity in terms 
of determinants. 
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A Equations of motion for spherically symmetric and asymptotically 
flat backgrounds 

Here we give details of the equations used in Section 2 to describe spherically symmetric solutions with 
the theory defined by (2.19). Using the Ansatz 



ds 2 = -N(r) 2 dt 2 + F{r)~ l dr 2 + r 2 H(r)~ 2 dQ 2 
in the action (2.19), we obtain the potential 
r 2 sm(9) 



7 gU = - Z~^ 2 y |2 \yf{{2H - 3)N + 1) + H 2 N + H(2 -6N) + 6N -3 

-6a 3 {H - 1) \y/F{(H - 3)N + 2) - 2HN + H + AN - 3 
-24a 4 (l - VF)(1 - H) 2 (l - AT) J. 
Varying it with respect to gives the equation of motion 

= 3a 3 mV(l - H)H 2 \Vf(H - 3) - 2H + a] - 12a 4 mV (y/f - l) (1 - H) 2 H 2 



r[F - 6m 2 r ) + 3m V VF + F 



+ rH rFH + 2F [rH + 3H 



-H< 

-5r 2 F (#J " + 2m V (VF - 3j H 3 + H 4 (m 2 r 2 + l) , 

where 'dot' denotes derivatives with respect to r. Varying with respect to F gives 

= -3a 3 mV(l - H)H 2 [H(2N — 1) — AN + 3] — 12a 4 mV(l - H) 2 H 2 (1 - N) 



r 2 FN (hY + H 2 



+2rFHH [rN + N 
+H^N (m 2 r 2 + l) + 2m 2 r 2 H 3 (l - 3N), 



N (6m V — F) — r (2FN + 3m 2 r 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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and finally by varying with respect to H, one gets 

^/F(H(2N - 1) - 3JV + 2) + H(2 - 3N) + AN - 3 
rFN + NF + 2F (rN + n) + 2m 2 r\ f F{3N - 1) - 12m 2 rN + 6m 2 r 



= Qa^m 2 rH 2 



+H rNFH + 2F rNH + rHN + 2NH 



)) 



4rFA [ HJ 



+2m 2 rH i 



VF-3)N + 1 + 24a 4 m 2 r (Vf - 1 ) (H - 1)H 2 (N - 1 



(A.5) 



Instead of using the last equation, it is simpler to use a constraint based on the (second) Bianchi identity. 
This can be achieved by taking a combination of the previous field equations which leads to \/^G^ v = 
for the Einstein piece, where G„v is the Einstein tensor. The constraint constructed in this way is 







1 



rHN 



H (N(4rH + 6) - 2rH + 3WV - 4) + 2r(2 - 3N)H 



-3a 3 



+rH 3 N + H 2 (-4rN - AN + 2) + 2H 2 (H(2N - 1) - 3N + 2) 



+12a 4 (l - H) 



2r(N - l)H + rH 2 N - H(rN + 2N-2) + 2H 2 (N - 1) 



+VF 



-H 



2N(rH + 3) + 3rN - 2 + 2(3iV - l)rH + 2H 2 (rN + N) 



-2H 2 [(H -3)N + 1] 

We would like to study perturbations about flat space, hence 

N = l + n, F = l + f, H = l + h, 



(A.6) 



(A.7) 



and study linear perturbations. However, n and / are small, like in GR, but h could be large since this 
is where the scalar graviton has an influence. Therefore, we will keep higher orders in h and truncate the 
equations to first order in n and /. It is more convenient to introduce a new radial coordinate p = r/H, 
so that the linearised metric is expressed as 



ds 2 



[l + 2n)dt 2 + (1 - f)dp 2 + p 2 dn 2 . 



(Ai 



The change of coordinates fixes / in terms of /, as f(p) = f(r(p)){d p [pH(p)]} 2 , while we have the 
freedom to choose n and h in terms of n and h. For simplicity, we pick n{p) = n(r(p)) and h(p) = h(r(p)). 
Therefore, we can drop the tildes of n and h from now on, and / simplifies to 



l-f(p) = [l-f(r(p))}[l + h(p)+phY, 



(A.9) 



where a prime denotes a derivative with respect to p. Then, equation (A. 3), in the new variable p and 
to leading order in f(p) but keeping all orders in h(p), reads 

= 2/ + 2pf + m 2 p 2 ( [1 - 2(3a 3 + l)h + 3(a 3 + Aa 4 )h 2 } [(2 + f)pti + (1 + h)f 

+6h[l - (3a 3 + l)h + (a 3 + 4a 4 )/i 2 ]) . (A.10) 

Equivalently for equation (A. 4), one gets 

= -/ - 2pn' + m 2 p 2 (n - 2[1 + n + (3a 3 + l)n]h + [(3a 3 + l)(n + 1) + 3(a 3 + 4a 4 )n]h 2 ) . (A.ll) 
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Finally, for the constraint equation (A. 6), the expression to first order in / and n, but to all orders in h 
is 

= pri [-1 + 2(3a 3 + l)h - 3(a 3 + 4a 4 )/i 2 ] - f[l - (3a 3 + l)h]. (A.12) 
These last three equations are the ones used in Section 2 to derived the Vainshtein mechanism. 



B General exact solution 

In order to construct exact solutions to the Lagrangian (2.19), we use the following Ansatz 

ds 2 = -C(r) dt 2 + 2D(r) dtdr + A(r) dr 2 + B(r) dfl 2 , (B.l) 

with 

B(r) = b r 2 , 

C(r) = c + ^ + c 2 r 2 , 

A{r) + C{r) = Q , 

D 2 (r) + A{r)C(r) = A . (B.2) 
The constant parameters bo, Co, c\, C2, Qo, Aq, can then be fixed using Einstein equations 

G^p = 8irG T^y (B.3) 

with T^v = — -^=^g-, and Ck given in Eq. (2.18). The combination C(r)G rr + A(r)G t t = fixes 
uniquely bo to be 

'1 + 6a 3 + 12a 4 + T ± X 2 



3(1 + 3q 3 +4a 4 ) 



(B.4) 



where 



r± = ±Jl + 3a 3 + 9a§ - 12a 4 . (B.5) 
By requiring that the 1/r 3 term in the Gtt equation vanishes, we obtain the condition 

* = fr- (B ' 6) 

and using the rest of Einstein's equations leads to a solution for the remaining coefficients. These are 
Anm 2 



c-2 



96 (l + 3a 3 + 4a 4 ) 2 L 
Qo ' 



1 - 2r± + 4a 4 (2r ± - 7) + 2a 3 (l - 18a 4 - 2T±) + a|(15 - 6T±) + 18a§ 



816o(l + 3a 3 + 4a 4 ) 4 



8 + 8r ± + 81A 12 + a 3 [10 + 81A + 2 6 3 3 al(2 + 3A ) + 9r± 

+2160^(1 + 18A + 4T ± ) + 6a 4 (17 + 162A + 22r ± )]r ± 
+27a| [27 + 162A + 288a^(5 + 9A ) + 20r ± + 8a 4 (29 + 162A + 26r±)] 
+2 7 3 3 a^(-l + 6A + 2T±) + 2 8 3 4 a^(l + A ) + 81a|(41 + 81A ) 
+54al(41 + 162A + 24a 4 (14 + 27A ) + 20r ± )r ± 

+144al(l + 54A + 8r±)r± + 48a 4 (2 + 27A + 3r±)l . (B.7) 

Notice that there are two branches, depending on the sign choice in bo (see equations (B.4) and (B.5)). 
For a 3 = — 4a 4 , the upper (lower) branch solution must be taken for a 4 < 1/12 (a 4 > 1/12). The solution 
for Qo exists only if Qo satisfies the condition 4VAo + 2Qo > 0. Due to this condition, for q 3 = — 4a 4 , 
the solution is valid only for a 4 in the ranges a 4 < 1/12 and a 4 > 1/8. 
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